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In the past years a good comprehension of the infrared gluon propagator has been achieved, with 
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a large lattice spacing. The interplay between volume effects and lattice spacing effects has not 
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propagator in the infrared region. Furthermore, we provide infinite volume extrapolations which 
take into account the finite volume and finite lattice spacing. 

We also report on preliminary results for the gluon propagator at finite temperature. 
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1. Introduction and motivation 

In a Quantum Field Theory, the knowlegde of all Green's functions allows a complete de- 
scription of the theory. In QCD, propagators of the fundamental fields encode information about 
non-perturbative phenomena, like confinement and dynamical chiral symmetry breaking. 

Here we study the gluon propagator in Landau gauge, at both zero and finite temperature, 
using lattice simulations. 

2. The gluon propagator at zero temperature 

In this section we study the gluon propagator using lattice simulations for various lattice vol- 
umes and lattice spacings. Although there is some support for D(0)=0 [§, ||], recent large volume 
lattice simulations, close to 20 fm, have claimed a finite non-vanishing gluon propagator at zero 
momentum In order to be able to simulate such large volumes, the reported simulations used 

a large lattice spacing ~ 0.2 fm. Since there is no systematic study about the effect of such a large 
lattice spacing in the propagator, our first goal is to investigate the interplay between volume effects 
and lattice spacing effects. Furthermore, we also consider the extrapolation of the propagator to the 
infinite volume limit. 

2.1 Lattice effects 

In figure [j], we show our results for the gluon propagator, renormalized at [i = 4GeV — see 
[f]] for details and the lattice setup 1 . For comparison purposes, we also show data obtained by the 
Berlin-Moscow- Adelaide collaboration [§]. Note that data in the same plot has the same lattice 
spacing, while varying the lattice volume. The plots show that, in the infrared region, the gluon 
propagator decreases as the lattice volume increases. 

Whereas in figure || we compare data with the same lattice spacing, in figure || we plot data 
with similar physical volumes. This allow us to study how the propagator changes with the different 
lattice spacings, keeping a constant physical volume. We have considered 4 different volumes, 
namely ~3.3,~4.6,~6.6, and ~8.1 fm. 

The first thing to note in figure |2] is that, for momenta above ~900 MeV, the propagator is 
well-defined, in the sense that the data define a unique curve. We can therefore claim that the 
renormalization procedure has been able to remove lattice artifacts for the high momenta region. 

Furthermore, figure ^| shows that the simulations performed with smaller j3 values (i.e. larger 
lattice spacings) underestimate the gluon propagator in the infrared region. Indeed, by comparing 
figures |] and || one can conclude that the corrections due to the finite lattice spacing seem to be 
larger than the corrections due to the finite volume. Moreover, from the above observations one 
can claim that the data from the Berlin-Moscow-Adelaide collaboration provide a lower bound for 
the gluon propagator in the continuum. 

'Simulations in this section have been performed with MILC code fm. 
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Figure 1: Renormalized gluon propagator for ji = 4 GeV for all lattice simulations. 

2.2 Zero momentum gluon propagator 

Figure || shows the zero momentum gluon propagator as a function of 1 /L. For each set of 
D(0) values with the same j8, one can consider the extrapolation to the infinite volume. Here we 
consider the ansatz 

D(0)= C -+D oo (0) (2.1) 

The results for D„(0) are 8.43±0.61GeV~ 2 (x 2 /d.o.f. = 2.6) for j8 = 5.7, 8.79±0.64GeV- 2 
(X 2 /d.o.f. = 1.7) for j8 = 6.0, and 9.72 ± 0.25GeV~ 2 {% 2 /d.o.f. = 0.1) for j3 = 6.2. In what 
concerns the data from the Berlin- Adelaide-Moscow collaboration, the fit provides Doo(0) = 6.1 ± 
1 .4GeV~ 2 (x 2 /d.o.f. = 1 .3). Note that this linear extrapolation does not provide a coherent picture 
of all the data sets. However, one can claim a Dco(0) in the range 6-10 GeV -2 . Note, however, 
that the results reported in [Q, ||] do not exclude completely a vanishing gluon propagator at zero 
momentum. 

2.3 Extrapolation to the infinite volume limit 

Having in mind the extrapolation of our results to the infinite volume limit, we consider fits of 
the lattice data to the functional form 



= Z . ^ 2 |f' . . (2.2) 
q 4 +M 2 q 2 +M* 



Note that for Z=l, the above expression is the tree level prediction of the so-called refined 
Gribov-Zwanziger action and, as shown in [[7|], it describes the infrared lattice gluon propagator up 
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Figure 2: Comparing the renormalized gluon propagator at fi = 4 GeV for various lattice spacings and 
similar physical volumes. 



j3 


Z 


M\ 


Ml 


M\ 


D(0) 


5.7 


0.821(10) 


4.09(17) 


0.558(36) 


0.380(11) 


8.84(45) 


6.0 


0.830(13) 


4.01(16) 


0.565(46) 






6.2 


0.83333(17) 


4.473(21) 


0.704(29) 


0.3959(54) 


9.42(14) 



Table 1: Extrapolation to the infinite volume limit. D(Q) is given in GeV 



to momenta ~ 1.5GeV. The extra parameter Z allows to fit the lattice data from up to 4GeV — 
see [§] for details. Then we combine all volumes for a given /3 and perform a linear extrapolation 
in 1/L to the infinite volume of each parameter independently. The extrapolated parameters can 
be seen in table [l[ All extrapolations have a % 2 /d.o.f. below 1.25 with the exception of M\ for 
j3 = 6.0; for this reason, no information is given about this parameter in table |] for the j3 = 6.0 
case. In figure Q we show the extrapolated propagators and compare with the largest lattice volume 
available in each case. 



3. The gluon propagator at finite temperature 

In this section, we focus on the calculation of the gluon propagator at finite temperature. On 
the lattice, finite temperature is introduced by reducing the temporal extent of the lattice, i.e. we 
work with lattices lf s x L t , with L t <C L s . The temperature is defined by T = 1 jah,. 

At finite temperature, the Landau gauge gluon propagator is described by two tensor structures, 
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Figure 3: Linear extrapolation of D(0) to the infinite volume. 
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Figure 4: Extrapolated propagators to the infinite volume. 
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Temp. (MeV) 




L s 


U 


a [fm] 


1/a (GeV) 


121 


6.0000 


64 


16 


0.1016 


1.9426 


162 


6.0000 


64 


12 


0.1016 


1.9426 


243 


6.0000 


64 


8 


0.1016 


1.9426 


260 


6.0347 


68 


8 


0.09502 


2.0767 


265 


5.8876 


52 


6 


0.1243 


1.5881 


275 


6.0684 


72 


8 


0.08974 


2.1989 


285 


5.9266 


56 


6 


0.1154 


1.7103 


290 


6.1009 


76 


8 


0.08502 


2.3211 


305 


5.9640 


60 


6 


0.1077 


1.8324 


305 


6.1326 


80 


8 


0.08077 


2.4432 


324 


6.0000 


64 


6 


0.1016 


1.9426 


486 


6.0000 


64 


4 


0.1016 


1.9426 



Table 2: Lattice setup used for the computation of the gluon propagator at finite temperature. 

Df v (q) = 8 ab (p£ v D T (q 4 ,q) +P^ v D L (q 4 ,q)^ (3.1) 
where the transverse and longitudinal projectors are defined by 

^v = (1-5 M 4)(1-M(V-^) - PU=(^- qj f)- p U (3-2) 
the transverse Dj and longitudinal Dl propagators are given by 

Dr{q) = 2v{N l 2_ l) ((Af (?)Af (-?)) - |f <A| (q)A a 4 (-q))^ (3.3) 



Dl{q) = V{n]-\) ( 1 + |( A 4(?M2(-?))J (3-4) 

In table || we show the temperatures we have simulated up to now 2 , for a fixed physical spatial 
volume ~ (6.5fm) 3 . For the determination of the lattice spacing we fit the string tension data in 



[10], using the functional form used in [11], in order to have a function a(j3). 

We resume the results obtained up to date in the 3d plots shown in figure ||. We see that 
the transverse propagator, in the infrared region, decreases with the temperature. Moreover, this 
component exhibits a turnover for small momenta. The longitudinal component increases for tem- 
peratures below T c ~ 270 MeV. Then the data exhibits a discontinuity around T c , and the propagator 
decreases for T >T r . 



Simulations in this section have been performed with the help of Chroma library [g|; the FFT transforms have been 
done with the PFFT library |p|. 
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Figure 5: Longitudinal (left) and transverse (right) gluon propagator as a function of momentum and tem- 
perature for a ~ (6.5fm) 3 spatial lattice volume. 
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